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La presents invention concerne un procede pour 
determiner les effets d'une sollicitation mecanique sur 
un objet const itue d'un materiau solide. 

Dans de nombreux domaines techniques, on a besoin de 
determiner 1 1 incidence de sollicitations mecaniques sur 
des* pieces, notamment des pieces metalliques, soit afin 
d'optimiser des procedes de fabrication comportant des 
operations de mise en forme, soit afin de prevoir le 
comportement en service de pieces dans des conditions 
complexes • 

Les procedes de mise en forme conceraes sont les 
procedes de mise en forme par deformation plastique tels 
que 1 1 emboutissage de toles, I'etirage de tubes, le 
f orgeage de pieces de mecanique . Ces procedes peuvent 
avoir une incidence importante sur les proprietes du 
materiau const itutif des pieces concernees, et il est 
important de pouvoir prevoir cette incidence de fagon^S. 
permettre d'evaluer a priori le comportement en service 
des pieces obtenues . 

Les sollicitations en service des pieces sont 
di verses et correspondent a toutes les sollicitations que 
1 1 on peut rencontrer, par exemple a l r effet d 1 efforts 
permanents, l f effet d 1 efforts appliques progressivement a 
des vitesses plus ou moins elevees, d 1 efforts appliques 
de fa<?on periodique, d 1 efforts appliques dans des 
conditions de temperature ambiante ou bien d 1 efforts 
appliques dans des conditions de temperature soit elevee 
soit trds basse. 

Afin d'evaluer les effets de ces sollicitations 
mecaniques sur des pieces, on utilise de fac?on connue des 
methodes de calcul numerique, en mode elastique ou 
anelastique. Ces methodes de calcul sont, notamment, les 
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methodes par differences finies et les methodes de type 
Galerkin reposant sur une decomposition sur des fonctions 
de base dont un type particulier est la methode par 
elements finis. Par "mode anelastique" , on entend tout 
5 mode de comportement qui n'est pas elastique reversible, 
c'est a dire tout mode de comportement plastique, 
plastique dependant du temps, ou viscoplastique, 
plastique dependant de la temperature, plastique 
resultant de transformation de phases. Ces calculs 
10 utilisent une modelisation du comportement du materiau 
sous I'effet des contraintes, representee par une loi 
reliant les contraintes aux deformations. De fagon 
classique, on utilise des lois de comportement 
macroscopiques qui correspondent a des representations 
15 phenomenologiques des phenomenes de deformation, c 1 est-a- 
dire des representations qui ne tiennent pas compte des 
phenomenes physiques f ondamentaux sous- j acents . De tels 
modeles ont I'avantage d'utiliser un nombre de variables 
relativement limite ce qui conduit a des calculs 
20 suffisamment peu volumineux pour pouvoir Stre utilises 
dans les processus de conception industriels. En 
revanche, ce type de modele presente 1 1 inconvenient de ne 
pas permettre de modeliser convenablement 1 1 ef f et de la 
microstructure des materiaux sur leurs proprietes, et en 
25 particulier sur les proprietes anisotropes, et plus 
encore sur devolution de ces proprietes anisotropes sous 
I'effet des sollicitations et des deformations. 

Afin de rem<§dier a cet inconvenient, il a ete propose 
d'utiliser des modules de comportement de materiaux, dits 
30 "modeles polycristallins " , dans lesquels le materiau est 
represents par une pluralite de grains ayant la structure 
cristallographique et les proprietes d'un grain du 
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materiau. Chaque grain est caracterise notamment par son 

orientation et par les systemes de glissement propres & 

i 

sa structure cristallographique . Dans ces conditions, le 
nombre de systemes de glissements est important puisqu'il 
va de douze pour un cristal . cubique a faces centrees a 
trente pour un cristal hexagonal compact . De tels modeles 
ont I'avantage de permettre une representation fine du 
comportement du materiau, mais 1 1 inconvenient de 
ntcessiter 1 utilisation d'un nombre tres important de 
variables. Ce grand nombre de variables conduit a des 
temps de calculs absolument prohibitifs qui ne permettent 
pas 1 1 utilisation de tels modeles dans un contexte 
industriel . Enfin aucun de ces modeles ne permet de 
prendre en compte les endommagements microscopiques des 
materiaux, et . done ces modeles ne permettent pas de 
prevoir les ruptures en service. r " 

Le but de la presente invention est de remedier & c^s 
inconvenients en proposant un procede pour determiner les 
effets d.'une sollicitation mecanique sur un objet 
const itue d'un materiau solide qui permette de simuler de 
fagon fine devolution des proprietes d'un 'materiau 
anisotrope, en particulier d'un materiau polycristallin, 
en utilisant un nombre suffisamment restreint de 
variables pour que les temps de calcul soient 
raisonnables et permettent de faire des simulations dans 
un contexte industriel, tout en permettant de prendre en 
compte les endommagements du materiau. 

A cet effet, 1 1 invention a pour objet un procede pour 
determiner les effets des sollicitations mecaniques sur 
une piece constitute d'un materiau solide, selon lequel 
on calcule les deformations et les contraintes engendrees 
par la sollicitation mecanique en une plural ite de points 
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de la piece au moyen d'une methode de calcul num6rique en 
mode anelastique, dans lequel le comport ement du materiau 
solide est represents par un modele de comport ement 
microscopique polycristallin utilisant une pluralite de 
briques "grain" ayant chacune une orientation et une 
fraction volumique, dont les deformations sont 
determinees a partir d'une pluralite de systemes de 
glissement propres au materiau solide, la trace du 
tenseur des deformations plastiques microscopiques de 
10 chaque brique "grain" <§tant nulle, les orientations 
distinctes et les fractions volumiques correspondant §. 
des briques "grain" etant representatives de 
1 1 anisotropie du materiau solide. Selon 1' invention, le 
nombre de briques "grains" est inferieur ou egal a dix et 
15 le nombre de systemes de glissement utilises dans le 
module est inferieur ou egal a six. En outre, le modele 
microscopique polycristallin peut eventue 11 ement utiliser 
jusqu'S. sept briques "joints de grains" ayant chacune une 
orientation et une fraction volumique, et dont les 
20 deformations sont determinees a partir d 1 au moins un 
systeme d ! ouverture, la trace du tenseur des deformations 
plastiques microscopiques de chaque brique "joint de 
grain" etant positive. Le module peut egalement utiliser 
eventue 11 ement une brique "trou" ayant une fraction 
25 volumique variable, dont les deformations plastiques 
sont purement volumiques et dont la trace du tenseur des 
deformations plastiques microscopiques est positive. 

De preference, les six systemes de glissement 
associes a une brique "grain" def inis dans un . repere 
30 orthonorme associe a la brique "grain" def inis par les 
normales au plan de glissement et par les directions de 
glissement sont tels que le double de la matrice 
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d 1 orientation associee a chague systeme de glissement 
n'est const itue que de "0", de "-1" et de "1" . 

Pour un materiau isotrope S structure cubique les 
deformations sont determinees de preference a partir de 
5 deux families de trois systemes de glissement, une 
premiere famille correspondant aux faces du cube 
cristallin, et la deuxieme famille correspondant aux 
plans de 45° de ces faces. De preference le nombre de 
briques "grain" est. egal a sept et les briques "grain" 
10 sont obtenues par symetrisation de deux briques "grain" 
primitives situees sur les grands cercles a 45° de la 
figure de poles. 

On peut en outre determiner 1 1 evolution de la texture 
du materiau en calculant la rotation Q du reseau 

cristallin de chaque brique "grain" par rapport a un* 
repere corotationnel Lorsque le materiau est isotrope a 
structure cubique, le nombre de briques "grain" est alors' 
de preference egal a dix. 

On peut egalement determiner 1 ' avancement d'une 
transformation de phase.. Pour cela, on peut attribuer a 
chaque brique "grain" , des caracteristiques f onction de 
1 1 avancement de la transformation de phase. 

D'une fagon generale, les orientations et les 
fractions volumiques des briques "grain" peuvent etre 
25 determinees par ajustement a partir d ! essais de traction 
biaxiale et de traction torsion. 

En 1" absence de brique "trou", le modele de 
localisation pour le passage du niveau macroscopique au 
niveau microscopique pour une brique "grain" , permettant 
30 de calculer le tenseur des contraintes microscopiques cr g 

dans la brique "grain" en fonction du tenseur des 
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contraintes macroscopiques S et du tenseur des 
deformations plastiques micros copiques dans la brique 
"grain" et du tenseur E P des deformations plastiques 
macroscopiques s ! ecrit : 

5 

a =E+2^(1-B) [T|(l® l)+al] (E p .6 p ) 

Dans cette formule a et r) sont deux coefficients 
dependant d'une part du materiau et d'autre part de la 
valeur maximale atteinte au cours de la deformation par 
10 le deuxieme invariant de von Mises de la deformation 
plastique macroscopique . 

De preference, lorsque les calculs sont des calculs 
par elements finis, on associe une seule brique "grain" 
et event ue 1 lement une brique "trou" ou une brique "joint 
15 de grain" a chacun des points d 1 integration de la methode 
de calcul par elements finis ; la repartition des briques 
entre les points d ! integration de la methode de calcul 
etant effectuee de telle sorte que cette repartition soit 
homogene . 

20 La methode de calcul par elements finis utilisee peut 

§tre une methode de calcul me c an i que ou thermo-mecanique 
faisant intervenir le temps, 

Dans ce cas, lorsqu'on utilise une brique "trou", la 
fraction volumique de la brique "trou" n'est pas nulle a 
25 1 1 instant initial des calculs. 

Lorsque le nombre de briques "joint de grain" ou de 
briques "trou" est superieur ou egal a un, on peut 
definir un critere tel que, lorsqu'en un point de I'objet 
le critere est satisfait, on attribue au point considere 
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une resistance mecanique du materiau nulle, de fagon a 
simuler la presence d*une fissure. 

L'objet en materiau solide est par exemple une piece 
metallique soumise a une sollicitation mecanique associee 
a une operation de mise en forme par deformation 
plastique, notamment telle que 1 1 emboutissage, l'£tirage, 
le laminage, le cintrage ou le f orgeage . Cette liste de 
procedes de mise en forme par deformation plastique n'est 
pas limitative. 

L'objet peut egalement etre une piece metallique 
soumise a une sollicitation mecanique correspondant * a 
1 1 utilisation de la pi§ce. L 'utilisation de la piece 
pouvant engendrer notamment des phenoraenes de fluage, de 
fatigue, ou de f atigue-f luage . 

Le procede selon 1' invention peut etre utilise pour" 
evaluer l r effet des procedes de mise en forme, tels que 
1 1 emboutissage de toles ou 1' extrusion d T un tube, sur les 
proprietes des. materiaux. Les resultats obtenus peuvent 
etre utilises par exemple dans des modeles de simulation* 
de tests de resistance aux . chocs de vehicules 
automobiles. 

Le procede selon 1' invention peut egalement etre 
utilise pour simuler la fabrication de boites metalliques 
telles que des canettes, obtenues par emboutissage, de 
fagon a en, determiner les proprietes d'emploi en tenant 
compte des phenomenes tels que les cornes d 1 emboutissage 
resultant de 1 1 anisotropie du materiau. 

Le procede selon 1 ' invention peut aussi etre utilise 
.pour simuler la fabrication de tubes d'echangeurs de 
chaleur en alliage de nickel pour les generateurs de 
vapeur ou pour simuler la fabrication de gaines en 
alliage de zirconium pour combustible de centrales 
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L 1 introduction de briques "joint de grain" et de 
briques "trou" dans le modele permet de simuler les 
endommagements microscopiques du materiau et done de 
prevoir les ruptures ou les risques de rupture en 
service- De ce fait, le procede selon 1 1 invention peut 
<§galement §tre utilise pour etudier le comportement et 
evaluer la duree de vie des installations. En particulier 
dans le domaine des centrales nucleaires pour prevoir la 
duree de vie de composants tels que les cuves, les 
generateurs de vapeur, les tuyauteries, les gaines de 
crayons combustibles, tous composants soumis a des 
charge s comp 1 exe s . 

Pour les centrales thermiques classiques, le procede 
peut etre utilise pour prevoir la duree de vie 
d 1 equipements soumis a des phenomenes de f luage . 

Dans 1 1 aeronaut ique, le procede peut etre utilise 
pour la prevision de la duree de vie des disques et des 
ailettes des reacteurs soumis a des phenomenes de fluage 
ou de f at igue- fluage . Le fluage est le mode 
d 1 endommagement d'une piece soumise a un chargement 
maintenu pendant une longue duree. La fatigue est le mode 
d 1 endommagement d'une piece soumise a un grand nombre de 
cycles de chargement successifs. La fatigue -fluage est la 
combinaison de ces deux modes d 1 endommagement , par 
exemple lorsque les cycles de chargement comportent un 
maintien de longue duree a contrainte 61evee . 

Ces utilisations ne sont donnees qu'a titre d 1 exemple 
et l ! homme du metier saura adapter le procede S toute 
situation particuliere faisant intervenir des phenomenes 
de plasticite ou de fluage. 

L l invention va etre maintenant decrite plus en detail 
en regard des figures annexees dans lesquelles : 
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la figure 1 represente les six systemes de 
glissement utilises par le modele polycristallin 
selon 1 1 invention. 

la figure 2 represente les six systemes de 
glissement s utilises pour representer le 
comportement d'un cristal cubique . 

- la figure 3 est la figure de poles des briques 
"grain" representant un materiau isotrope. 

la figure 4 represente les six systemes de 
glissement utilises pour representer le 
comportement d'un cristal hexagonal compact. 

- la figure 5 est la figure de poles des briques 
" j oint de grain" d 1 un modele d 1 endommagement 
intergranulaire . ■ • * 

- la figure 6 represente 1 1 affectation des briquesy 
aux points de Gauss d'un element de calcul par..: 
elements finis. 

- la figure 7 represente 1 ' affectation des briques,-; 
aux points de Gauss d f un element d 1 une coque. 

- la figure 8 est la figure de poles des briques 
"grain" utilisees pour representer la texture d'un 
materiau initialement isotrope. 

- la figure 9 montre les resultats de calculs 
correspondant a un essai de traction biaxiale sur 
un tube en Zircaloy recristallise . 

Alors que dans ■■ les modeles de comportement 
polycristallins connus on utilise autant de systemes de 
glissement qu ! il en existe dans le cristal reel, de douze 
a trente selon la nature du cristal, les inventeurs ont 
decouvert de facpon nouvelle qu ! urie bonne representation 
de la deformation d'un cristal pouvait etre obtenue en 
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utilisant seulement six systemes de glissement distincts 
bien choisis. 

Alors que les modeles polycristallins connus 
considerent qu'il est necessaire d'utiliser autant de 
5 grains qu'il y a d 1 orientations effectives des grains 
dans le materiau reel, les inventeurs ont decouvert de 
fagon nouvelle, qu'il est possible d ! obtenir une bonne 
representation du comport ement du materiau en utilisant 
un nombre de grains restreint correspondant a au plus dix 
10 orientations differentes des grains. 

Les inventeurs ont egalement decouvert de fagon 
nouvelle qu'il etait possible dans tous les cas de 
chargement , notamment cycliques , et pour les briques 
"grains" comme pour les briques "joint de grain", 
15 d'obtenir une bonne representation des materiaux en 
utilisant des modeles de changement d'echelle permettant 
de faire le passage de 1'echelle macroscopique a 
1'echelle microscopique et inversement qui ne fasse pas 
intervenir de tenseur supplement aire mais qui, a la place 
20 de ces tenseurs, utilise des parametres scalaires. 

En outre, les inventeurs ont constate de fagon 
nouvelle qu'il n 1 etait pas indispensable d ! affecter 
1 1 ensemble des grains du modele de comportement a 
1 1 ensemble des points de calcul, et dans le cas du calcul 
25 par la methode des elements finis, a 1' ensemble des 
points d 1 integration de la methode par elements finis, 
mais qu'il etait possible d'obtenir une precision 
excellente en affectant un nombre r^duit de grains & 
chacun des points de calculs, et dans le cas de calcul 
30 par methode des elements finis, un grain a chacun des 
points de calcul de la methode consideree. 
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Enfin, les inventeurs ont decouvert de fagon 
nouvelle, qu'il etait possible de simuler 1 1 endommagement 
a l'echelle microscopique en utilisant des modeles de 
joint de grains simples qui permettent de simuler les 
phenomdnes d 1 endommagement intergranulaires ou d f utiliser 
des modeles de cavites pour modeliser 1 1 endommagement 
ductile des materiaux.' 

Ainsi, 1" invention se -distingue des methodes de 
calcul numerique connues par la nature du model e de 
comportement polycristallin de materiau, et par la fagon 
de 1 1 incorporer dans la methode de calcul numerique . 

Afin de bien faire comprendre 1 1 invention, on va 
considerer lecas particulier de la methode de calcul par 
elements finis. On va tout d'abord rappeler quelques 
notions sur les techniques de calculs par cette methode 
et sur les modeles de comportement de materiaux de typje 
polycristallin, puis on va expliquer comment, selon 
1' invention, on construit un modele de comportement 
polycristallin performant, et on va indiquer comment on 
20 incorpore ce modele de comportement polycristallin dans 
la methode de calcul par elements finis. 

De fagon classique, pour simuler par un calcul par 
elements finis le comportement d'un objet soumis a des 
sollicitations mecaniques, on construit une 

25 representation de 1 ! objet materiel considere a 1 1 aide 
d ! un maillage, et on identifie les sollicitations 
mecaniques exterieures auxquelles ' il va etre soumis. Puis 
de fagon connue en "elle-melme, on calcule par iteration, 
les contraintes et les deformations a des instants 
30 successifs separes par un pas de temps At. Ces calculs 
font intervenir un modele de comportement du materiau. 
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Dans le cas d l un calcul par elements finis utilisant 
un model e de comportement du mater iau du type "modele 
polycristallin" , a un instant determine t, on connait 
pour chacune des mailles du maillage, le tenseur des 
5 contraintes macroscopiques Z et le tenseur des 

deformations macroscopiques E qu'on peut decomposer en un 

tenseur des deformations elastiques E e et un tenseur des 

deformations plastiques B P et pour chacun des grains du 
modele polycristallin, le tenseur des deformations 
10 plastiques microscopiques . 

A 1 ■ aide d'un modele de localisation permettant le 
• passage de l'echelle macroscopique & l'echelle 

microscopique, on determine le tenseur des contraintes a 

=s 

dans chacun des grains du modele polycristallin. 
15 A l'aide d'un moddle de comportement anelastique tel 

qu'un modele plastique ou viscoplastique (de tels modeles 
sont connus en eux-memes de I'homme du metier), et des 
systemes de glissement du grain, on determine le tenseur 

.p 

des vitesses de deformation microscopiques de chaque 

20 grain. 

En integrant les vitesses de deformations sur un 
interval le de temps At, on calcule les deformations 
plastiques microscopiques & F a 1 1 instant t + At . 

=3 

En utilisant un module d 1 homogeneisation permettant 
25 le passage de l'echelle microscopique a l'echelle 
macroscopique, on determine alors les tenseurs des 
deformations E et des contraintes Z macroscopiques a 

1 1 instant t + At . On peut alors recommencer le calcul et 
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ainsi determiner Involution au cours du temps des 
contraintes et des deformations dans 1'objet <§tudie. 

Pour mettre en ceuvre une methode de calcul numerique 
du type de celle qui vient d'etre decrite, on utilise 
5 trois repdres: 

un repere global, ou repere du laboratoire, (I*, L 2/ 
L 3 ) qui est absolu et qui sert de reference pour la 
description de 1'objet a etudier. Selon la nature de ■ 
1'objet etudie, ce repere permet de reperer chaque point 
10 de 1'objet par ses coordonnees cartesiennes ou par ses 
coordonnees cylindriques , ou par ses coordonnees 
spheriques . 

des r.eperes locaux macroscopiques orthonormes, 
associ£s a chacun des points de I'objet etudie. Ainsi, a 
15 chaque point, et en particulier, a chaque point pour 
lequel on effectue un calcul, on associe un repere local 
. macroscopique. Ce repere est un repere orthonorme dont la 
position et 1 ' orientation par rapport au repere global 
peut evoluer au fur et a mesure que I'objet etudie se 
20 de forme. 

des reperes locaux microscopiques ortbionormes 
associes a chacune des "briques" du modele polycristallin 
et qui servent a decrire les deformations microscopiques 
de chaque brique . L 1 orientation du repere local 
25 microscopique d'une !, brique n associee a un point de 
I'objet etudie reprdsente 1 ' orientation de ce grain par 
rapport au repere local macroscopique du point considere. 

Par la suite on appellera (X X/ X 2 , X 3 ) un repdre local 
macroscopique et (xx, x 2 , x 3 ) un repere local 
30 microscopique. 

On va maintenant decrire le modele polycristallin 
utilise dans 1 ! invention. 
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Selon 1' invention, pour construire le modele 
polycristallin, on utilise des briques elementaires de 
trois natures differentes : 

des briques "grains" qui representent le 
comportement de grains cristallins ; 

~ les briques "joints de grains" qui representent le 
comportement de joints de grains ; 

- les briques "trou" qui representent des cavites. 
Les briques "grains" permettent de represent er le 
materiau en 1 1 absence d 1 endommagement microscopique . Les 
briques "joints de grains" permettent de representer des 
endommagements du type intergranulaire . Les briques 
"trou" permettent de representer les endommagements du 
type ductile. 

On va maintenant decrire d'abord les briques 
"grains", puis les briques "joints de grains" et enfin 
les briques "trous" . 

De fagon connue, a chaque brique "grain" qui est 
1' equivalent des "grains" utilises dans les modeles 
polycristallins connus, on associe des systemes de 
glissements, des lois de comportement du materiau pour 
chacun des systemes de glissement, une orientation de la 
brique "grain" par rapport a un repere lie a 1'objet que 
l ! on veut etudier, et enfin une fraction volumique f . 

Selon 1' invention, et quelle que soit la nature 
cristallographique du materiau solide conaidere , on 
attribue a la brique "grain" six systemes de glissement. 
De fagon connue, chaque systeme est d^fini par la normale 
a un plan de glissement et par une direction de 
glissement, ce qui permet d'associer S chaque systeme de 
glissement une matrice d 1 orientation . Dans un repere 
local microscopique orthonorme (x X/ x 2 , x 3 ) lie a la 
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brique "grain", selon 1 1 invention, on choisit les 
systemes de glissement de telle sorte que les doubles des 
matrices d 1 orientation, 7^ des systemes de glissement ne 

soient constitues que de 0, 1, -1. 

Si n s est le vecteur normal au plan du systdme de 
glissement d ! indice s, et Is 1© vecteur de la direction 
de glissement de ce systeme de glissement, on a par 
definition : 

. Avec ces definitions, les systemes de glissement, 
representes §. la figure 1, sont les suivants : 



n s ~ 


(i . 


0 


0) 


Is = 


(0 


1 


0) 
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ffis = 
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0 
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s = 4 n s = 


(1/^2 


1/V2 


0) 


Is = 


( -1/ Jo 


1/J2 


0) 




-i 


0 


0 


2 m s ~ 


0 


1 


0 




A 
U 


u 


u 


s = 5 n s = 


(0 


1/^2 


1/^2) 


Is = 


(0 


-1/-J2 

f v 


1/^/2 ) 




0 


0 


0 


2 m s = 


0 


-1 


0 




0 


0 
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s = 6 n s = 




0 


1/V2) 


Is = 


(1/V2 


0 


-1/^) 
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0 
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2 m s - 
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0 




0 


0 


-1 



A la figure 1, on a represents un cube ABCDEFGH dans 
le repere local microscopique (x 1; x 2/ x 3 ) . 

Le systeme de glissement s=l, correspond a un 
5 glissement selon le vecteur li, parallele a 1 1 axe x 2/ dans 
un plan tel que le plan CDEH qui est perpendiculaire au 
vecteur rii et parallele au plan (x 2/ x 3 ) . 

Les systemes de glissement s=2 et s=3 , sont 
comparables au systeme s=l, et correspondent a des 
10 glissement s paralleles aux axes x 3 et x 1; le long de plans 
paralleles aux plans (x 3 , x x ) et (x 1; x 2 ) , respect ivement . 

Le systeme de glissement s=4 correspond a un 
glissement selon le vecteur I4 paralldle a la diagonale DF 
du carre ADEF situe dans le plan (x x , x 2 ) , le glissement 
15 se faisant dans un plan tel que le plan DFGC 
perpendiculaire au vecteur n 4 . Ce plan est parallele a 
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l'axe x 3 et forme un angle de 45° avec les plans (x 1# x 3 ) 
et (x 2 , x 3 ) . 

Les systemes de glissement s=5 et s=6 sont 
comparables au systeme s=4 et correspondent & des 
glissements paralleles aux diagonales FB et BD, le long 
des plans FBCE et BDEG, respectivement . Ces plans sont 
paralleles aux axes x x et x 2 respectivement, et forment 
des angles de 45° avec les plans (x 2 , Xi) et (x 3/ x x ) pour, 
l'un, et avec les plans (x X/ x 2 ) et (x 3 , x 2 ) pour 1 1 autre . 

• 

Avec de tels systemes de glissement, si y 

represente la vitesse de glissement dans la direction l s 
du systeme de glissement s, le tenseur des vitesses de 
deformation plastiques de la brique "grain" s'ecrit : 



£ = 2 (m s y s )=± 



r. . .' 

• • • • 

7 7 7 ~7 , 



Pour determiner les vitesses de glissement y s sur 

chacun des systemes de glissement s, on utilise un modele 
de comportement qui utilise les relations de 
viscoplasticite classiques decrites par exemple dans La 
Revue Europeenne des elements, finis Vol 3 - n°4/l994, 
pages 515 a 541. Utilisation de modeles polycristallins 
pour le calcul par elements finis. Georges Cail let aud - 
Philippe Pilvin et qui font intervenir de fagon connue la 
scission resolue x B sur le systdme s. La scission r£solue 
x s s 1 ecrivant : 
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Dans cette equation, a est le tenseur des 

=s 

contraintes dans le grain considere. 

Pour modeliser le comportement d ! un materiau 
determine, il faut alors choisir des briques "grain" du 
modele polycristallin de telle sorte que ces briques 
representent d'une part le comportement d'un grain du 
materiau, compte-tenu de sa nature cristallographique 
reelle et d' autre part d'un agregat de tels grains, 
compte-tenu de la repartition reelle des orientations des 
grains dans le materiau, c'est-a-dire de la texture du 
materiau. 

Pour faire ce choix, on associe a chague systeme de 
glissement une scission critique telle que la brique 
"grain" obtenue ait un comportement identique a celui du 
cristal reel dont est constitue le materiau, puis on 
choisit le nombre de briques "grains", leurs orientations 
et leurs fractions volumiques de telle sorte que 
1' agregat obtenu ait un comportement bien representatif 
du matSriau reel . 

Pour representer 1 1 orientation des grains, on 
considere un repere local macroscopique (X x , X 2 , X 3 ) et un 
repere local microscopique (xi, x 2 , x 3 ) lie a chaque 
brique "grain". L 1 orientation d'une brique "grain" est 

definie comme etant 1 ! orientation de son repere local 
microscopique (x 1; x 2 , x 3 ) par rapport au repere local 

macroscopique (X X/ X 2/ X 3 ) . 

De fagon connue en elle-meme, on considere que le 

repere local microscopique (x X/ x 2/ x 3 ) se deduit du 

repere local macroscopique (X x , X 2 , X 3 ) par trois 

rotations : 
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- une rotation d 1 angle V[/ du repere (X a , X 2 , X 3 ) autour 
de l'axe X 3 qui conduit a un repere (X ! i, X' 2/ X 3 ) , 

- une rotation d 1 angle ({> autour de l'axe X'i qui 
conduit a un repere (X 1 !, X" 2 , X' 3 ) 

5 • - une rotation d' angle cp du repere (X'1, X" 2 , X' 3 ) 

autour de 1 1 axe X'3 qui conduit a un repere (X n x , X n, 2 , 
X' 3 ) ; ce dernier repSre est le repere local microscopique 
(x x , x 2 , x 3 ) . 

Le repdre local microscopique (xi, x 2 , x 3 ) est done 
10 defini par les trois angles de rotation *F, <j>, <p, qui ont 
permis de l'obtenir. Ces angles, sont habituellement 
appeles "angles d'Euler" . 

Cette methodologie de choix des briques "grains" ya 
etre illustree par deux exemples que I'homme du metier 
15 saura transposer a tout cas particulier. 

Le premier exemple, concerne le cas d ! un materiau 
isotrope const itue de cristaux de structure cubiqvie 
centree ou f aces-centrees tels que les aciers. 

Lorsque le materiau est constitue de cristaux de 
20 structure cubique, les ' systemes de glissements. peuvent 
etre regroupes en deux sous -families correspondant 
chacune a une scission critique. Ces deux families de 
systemes de glissement correspondent d'une part aux 
glissements sur les f aces ' du cube, et d 1 autre part aux 
25 glissements sur les plans a 45° de ces faces. 

Ces six systemes de glissement adaptes a un cristal 
cubique sont representes & la figure 2. Les systemes s=l, 
s=2 et s=3, correspondent aux glissements sur les faces 
du cube et correspondent a la meme scission critique, 
30 puisque, hormis 1 1 orientation, ces faces sont identiques. 
Les systemes s=4, s=5 et s=6, correspondent aux 
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glissements sur les plans a 45° par rapport aux trois 
directions des aretes du cube. Ces trois systemes 
correspondent a la meme scission critique, puisque, 
hormis 1 1 orientation, ces plans sont identiques, 

5 Afin de bien representer le caractere isotrope du 

materiau, on considere un modele constitue de sept 
briques grains dont les orientations sont choisies de 
fagon a couvrir au mieux une sphere. 

On utilise alors 7 briques "grain" choisies de telle 

10 sorte que la figure de poles correspondante soit 
symetrique . 

Pour cela, on construit un premier groupe de trois 
briques "grain" et un deuxieme groupe de quatre briques 
"grain" en procedant comme suit : 

15 - pour le premier groupe, on considere une premiere 

brique "grain" primitive d'indice g=l, obtenue par une 
rotation d 1 angle v F=ti:/4 autour de l ! axe X 3 . Pour obtenir 
un comportement identique selon les trois axes, on 
symetrise la figure de p61es en ajoutant deux briques 

20 "grain" d'indice g=2 et g=3 obtenues par rotation de 7t/4 
autour des axes X x et X 2 respect ivement . Les coordonn6es 
des poles de ces grains sont alors 0, ± 1, ± 1/V2 et 
± (a/2-1) . 

- pour le deuxieme groupe, on considere une deuxieme 
25 brique "grain" primitive dont les poles sont situes sur 
trois grands cercles passant par un des axes Xi et 
formant un angle de 45° avec les deux autres axes, et 
dont les coordonnees sont ± 0,2, ± 0,4 et ± 0,5. Puis on 
symetrise la figure de poles en ajoutant trois briques 
30 "grains" dont les poles se deduisent des poles de la 
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deuxieme brique "grain" primitive par symetrie par 
rapport aux axes Xi, X 2 et X 3 , 

On obtient ainsi sept briques "grains" auxquelles on 
attribue la fraction volumique f = 1/7 et dont les 
5 angles, ¥, <j>, et (p, exprimds en radians, sont : 

grains: 1 2 3 4 5 6 7 

n / 4 0 tc/2 0,464 2,034 .-2,678- 

4 0 tt/4 7t/4 0,841 1,231 '0,841 0,841 

9 0 0 0,464 -tc/4 -1,107 2,034 



Les sens de deux des trois axes xl, x2 et x3 peuvent 
etre inverses, et les trpis axes peuvent etre soumis & 
10 1 une permutation circulaire sans modifier les resultats du 
modele, Les angles du tableau correspondent a un choix 
particulier des sens des axes et de la permutation. 

La figure de p61es. obtenue est representee a la 
figure 3. Sur cette figure, les poles rep^res l x , I2/ I3 
15 representent les poles de la premiere brique "grain" . Les 
poles rep^res 2 Xl 2 2 , 2 3 representent les poles de la 
deuxieme brique "grain", et ainsi de suite. On voit en 
particulier que les pdles 4 3 , 5 2 , 61 et 7 3 ont pour 
coordonnees (+ 0,5 ; ' + 0,5). Les pSles 4 X/ 4 2 , 5 lt 5 3 , 6 2/ 
20 63, 7i et 7 2 ont pour coordonnees (+ 0,2 ou + 0,4) . 

Ce modele est universel pour les materiaux isotropes 
a structure cubique centree ou faces-centrees. 

A titre de deuxieme exemple, on considdre un materiau 
anisotrope ayant une structure hexagonale compacte. Un 
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tel materiau est par exemple du Zircaloy recristallise . 
Pour une telle structure cristallographique hexagonale, 
on considdre trois groupes de systemes de glissement 
correspondant a trois scissions critiques distinctes. Le 
premier groupe rassemble les systemes 2 et 3 qui 
correspondent au glissement basal, le deuxieme groupe 
rassemble les systemes 1 et 4 qui correspondent au 
glissement prismatique, et le troisieme groupe rassemble 
le systeme 5 et 6 qui correspondent aux glissements 
pyramidaux. 

Les systemes de glissement destines a representer un 
cristal hexagonal compact sont representes 3. la figure 4. 
Sur cette figure, les axes x X/ x 2 , x 3 constituent le 
repere orthogonal de reference de la brique "grain", 
c'est-l-dire le repere local microscopique . Le prisme a 
base hexagonale ABCDEFGHIJKL est une structure 
qui dans le cadre du modele a un comport ement identique a 
celui d'un cristal hexagonal, mais n*a pas la meme 
geometrie que celle d'un cristal de structure hexagonale. 
Les systemes de glissement : 

- s=2 repere P 2 sur la figure, et s«3, repere B 3 
correspondent au glissement basal ; il leur est 
attribue une premiere valeur de scission critique ; 

- s=l repere P x et s=4 repere P 4 , correspondent au 
glissement prismatique ; il leur est attribue une 
deuxieme valeur de scission critique ou deux valeurs 
distinctes . 

- s=5 repere tc 5 et s=6 repere % S r correspondent aux 
glissements pyramidaux ; il leur est attribue une 
valeur de la scission critique particuliere . 
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En outre, il est necessaire de choisir un certain 
nombre de briques "grains" ayant des orientations 
permettant de representer 1 ' anisotropie globale du 
materiau. Dans le deuxieme exemple, qui concerne en 
5 particulier un tube en Zircaloy recristallise, 
1 ! invention utilise six briques "grains" ayant chacune 
une fraction volumique de 1/6, dont les orientations par 
rapport au repere local macroscopique (Xi, X 2 , X 3 ) lie au 
point de l'objet considere, sont les suivantes : 



g = 
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2 . 
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71/ 2 
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7C-T 2 


<j> . 


tc/2 


71/ 2 
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Les orientations • de ces six grains dependent de 1 cinq 
parametres qui correspondent a cinq angles d'Euler *F X/ 
^2/ <Pi, <?2 et 93. Ces cinq parametres peuvent §tre 
15 determines par une methode d 1 identification numerique en 
comparant les resultats de calcul numerique. a des 
resultats d'essais de traction-torsion. Ces methodes 
d 1 identification numerique sont des methodes connues de 
l'homme du metier. 
20 Pour modeliser Involution de la texture d'un 

materiau polycristallin, c r est-a-dire l 1 evolution des 
orientations des grains, on considere de fagon connue la 
rotation Q du reseau cri stall in d ! un grain par rapport a 

un repere dit "corotationnel » . La vitesse de rotation du 
25 reseau cristallin est donnee par : 
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. A' 

ou apparaxt le tenseur de rotation : 

Le repere corotationnel suit la rotation moyenne de 
5 la matiere R . De fagon connue, sa vitesse de rotation par 
rapport au repdre du laboratoire (LI, L2 , L3) est 
calculee a partir des gradients des vitesses v des points 
materiels du corps solide d^formable : 

R:R T =( grad v~ grad T v)/2 

10 Dans I'exemple du materiau initialement isotrope 

constitue de 7 briques "grain", il peut etre interessant, 
pour certaines geometries de pieces et certains 
chargements particuliers comme le laminage d'une tole, de 
maintenir une symetrie orthotrope apres evolution de la 

15 texture. Pour cela et de fagon nouvelle, on dissocie les 
grains 1 a 3, obtenus par rotation de autour des 

trois axes X3 , XI et X2 du repere local macros cop i que, en 
six grains obtenus par des rotations de 7c/8 (grains 1 §. 
3) et 3tc/8 (grains 8 a 10) autour de ces axes. De plus, 

20 les axes x2 et x3 du repere a I'echelle microscopique du 
grain 8 sont inverses, les axes xl et x3 du grain 9 sont 
inverses et les axes xl et x2 du grain 10 sont inverses, 
ce qui est equivalent a des rotations d 1 angle n autour 
des axes xl, x2 et x3 respectivement pour les grains 8, 9 

25 et 10. On obtient ainsi dix briques "grains", dont les 
angles d'Euler ^ et (p et les fractions volumiques f 

sont : 
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grains 
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evolution de texture, est representee a la ■ figure 8. Pour 
obtenir une figure de poles orthotrope apres evolution de 

5 texture, les dix briques "grains" 1 a 10 peuvent etre 
completees par dix briques "grains" obtenues ' en 
substituant les unes aux autres les norraales aux plan's de 
glissement . n s et les directions de glissement [ s des 6 
systemes de glissement , ce qui revient a conserver les 

10 tenseurs d r orientation m , et §l inverser les tenseurs de 

rotation q . 

—s 

On va maintenant deer ire les briques "joint de 
grains", destinees a modeliser les endommagements 
intergranulaires lies a I'ouverture des joints de grains. 
15 Pour mod<§liser les endommagements intergranulaires , 

on utilise des briques "joint de grains" dont le 
comportement est represents par au moins un systeme 
d ' ouverture et e vent ue 1 1 ement au moins un systeme de 
glissement. 

20 Un systeme d 1 ouverture se distingue d'un systeme de 

glissement par le fait qu'il conduit a un tenseur de 
deformations plastiques dont la trace est positive. 
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Comme pour les briques "grains", les systemes 
d'ouverture cm de glissement sont definis par des 
matrices d x orientation m se deduisant des normales aux 

plans de deformation (glissement ou ouverture) « s et des 
5 directions de deformation [ s . 

A titre d'exemple, on considere les deux cas suivants 

1) la brique "joint de grain" comporte un systeme 
d 1 ouverture (s=l) et deux systdmes de glissement 
10 intergranulaire (s=2 et s=3) . Ces systemes sont definis 
de la facpon suivante : 



s = 
n s = 



(1 0 0) (1 0 0) 



(1 0 0) 



l s = (1 0 0) (0 1 0) (0 0 1) 
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2) La brique "joint de grain" comprend trois systemes 
15 d'ouverture (s=l, s=2 , s=3) et trois systemes de 
glissement intergranulaire (s=4, s=5, s=6) : 



s = 1 2 3 

n s = (1 0 0) (0 1 0) (0 0 1) 

1 6 - (1 0 0) (0 1 0) (0 0 1) 
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Le choix des briques "joint de grain" depend de la 
nature des sollicitations auxquelles est soumis l ! bbjet 
5 dont on veut etudier le comportement . 

Lorsque 1 1 obj et dont on veut etudier le comportement 
est soumis a des sollicitations dans une seule direction, 
XI , on peut considerer une seule brique "joint de grain" 
perpendiculaire a la direction XI, dont le . systeme 
10 d'ouverture correspond a la direction XI et les deux 
autres systemes de glissement correspondent a X2 et X3 . 

Pour un . chargement dans une direction quelconque , on 
peut soit. utiliser trois briques "joint de grain" 
perpendicul aires a chacune des directions XI, X2, X3 , 
15 soit utiliser une seule brique "joint de grain" 
comportant trois systemes d'ouverture et trois systemes 
de glissement. 

Enfin pour la description f ine - de 1 1 endommagement par 
une sollicitation multiaxiale, on peut considerer sept 
20 briques ' "joint de grain" ayant chacune un systeme 




d'ouverture et deux systemes de glissement , et une 
orientation definie par les deux angles d'Euler de la 
normale n du systeme d'ouverture. 

Les orientations des directions d'ouverture des 
5 briques "joint de grain" peuvent etre representees par 
une figure de poles reproduite a la figure 5, dans 
laquelle les sept orientations sont reper6es 1, 2, 3, 4, 
5, 6 et 7 . 

Les fractions volumiques des briques "joint de grain" 
10 sont choisies de telle sorte que leur fraction volumique 
totale soit egale & six fois le rapport de l'epaisseur 
mecanique du joint de grain a son diametre. A titre 
indicatif, cela represente quelques pour-cent. 

Pour decrire le comportement des joints de grains, on 
15 utilise des modeles visco-plastiques classiques qui 
permettent de determiner les vitesses d'ouverture des 
joints de grain. Ces modeles sont par exemple les modules 
decrits par Onck et Van der Giessen dans Journal of the 
Mechanics and Physics of Solids -Vol 47, 1999, pages 99- 
20 139- "Growth of an initially sharp crack by grain boundary 
cavitation" . 

On va maintenant decrire une brique "trou", destinee 
a modeliser 1 1 endommagement ductile. Une telle brique 
"trou" est caracterisee par une vitesse de deformation 
25 plastique purement volumique representee par les 
equations : 



1/3^ exp 



E p 
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p = exp(-3^) 

Em = f &m 

f = l-(l-f 0 )exp (-3EC) 
dans lesquelles : 

Dl e t a x sont des coefficients scalaires caracteristiques 
du materiau, 

5 s£ est la vitesse de deformation plastique volumique 
microscopique . 

2^ est la contrainte macroscopique moyenne, egale au 1/3 
de la trace du tenseur des contraintes . 

est la vitesse de deformation plastique macroscopique 
10 equivalente, au sens de von Mises. u 

E£ est la vitesse de deformation plastique macroscopique 
moyenne , 

f 0 est la fraction volumique initiale de la brique "trou" 

f est la fraction volumique de la brique "trou" a 
15 1 1 instant t . 

p est la masse volumique relative. 

Cette brique "trou" a une • fraction volumique f qui 

est une variable. Quand on ajoute une brique "trou" a un 

modele constitue de briques "grain" de fraction volumique 
20 f g chacune, on considere que la fraction volumique de 

chaque brique "grain" devient (1-f) f g . 

Les differentes briques "grain", "joint de grain" et 

"trou" qui viennent d'etre decrites, permettent de 
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modeliser le comportement microscopique du materiau. 
Mais, comme on l»a vu precedemment pour modeliser le 
comportement d'un objet, il est necessaire de pouvoir 
passer du niveau microscopique au niveau macroscopique et 
reciproquement . Pour cela, on utilise des modeles appeles 
"modeles de changement d'echelle" dont le principe 
general est connu de l'homme du metier, mais qui, dans le 
cadre de la presente invention, comportent des 
particularity qu l on va decrire maintenant. 

On considere trois situations selon que le modele ne 
comporte que des briques "grain 11 , que le modele comporte 
des briques "grain" et des briques "joints de grain" ou 
que le modele comporte des briques "grains" et une brique 
"trou" . 

Lorsque le modele ne comporte que des briques 
"grains", les relations permettant le passage de 
l 1 echelle macroscopique a l 1 echelle microscopique et 
inversement sont respect ivement : 

a) pour le passage de 1 ' echelle macroscopique a 
1 1 echelle microscopique : 

Zg" 2 + 2 V (l-p)a(E p -e P g ) 

avec : 

- = 1 + D(E P )™* 
a ~ 

et 

P " 15 l-v 

Dans ces equations , 

cr^ est le tenseur des contraintes microscopiques dans 
la brique "grain" d'indice g. 
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£ est le tenseur des contraintes macroscopiques . 

as 

est le tenseur des deformations plastiques 
macroscopiques . 

s P est le tenseur des deformations microscopiques 

dans la brique "grain" d'indice g. 

\x est le module de cisaillement du materiau. 
v est le coefficient de Poisson du materiau, 
D est un parametre propre au materiau, a just able a 

partir d'essais. 

(^l^est la valeur maximale du deuxi£me invariant de 

= eq 

von Mises du tenseur des deformations plastiques 
macroscopiques, atteinte anterieurement . 

b) Pour le passage de l'echelle microscopique a 
l'Schelle macroscopique : 

S= Z (fg£g) 




dans ces equations, 

f g est la fraction volumique de la brique "grain" 
d ' indice g, 

E* est le tenseur des deformations elastiques 
macroscopiques . 



32 



niveau macros copi que au niveau microscopigue et du 
passage du niveau microscopigue au niveau macroscopique 
sont modifiers de la facon suivante : 



a = £ + 2^(1 



-^77(l®l) + al 



(E p - s") 



avec 

i=i + D(^)r 

et 

77 — 



10 Enfin, lorsque le modele comporte des briques 

"grains" et une brique "trou", les relations de passage 
du niveau macroscopique au niveau microscopique et 
inversement sont les suivantes : 

a) pour le passage de l'echelle macroscopique a 

15 l'echelle microscopique : 
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a = 

=* 1-/ 



et 



2 + Dyf <r, exp 



Per, 



VA=1 = =. 



£ 0 = "(I - /) A a", exp 



Per. 



£0 est le tens eur des contraintes microscopiques dans 
la brique "trou". 

Di et o x sont des parametres propres au materiau. 

b) et pour le passage de l ! <§chelle microscopique a 
l'echelle macroscopique : 
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2 = fcr + <l-f) £ (fgcr ) 

E P =f6 P l +(l-f)' J (fg€ P ) 

5 Les briques telles qu'elles ont et<§ definies plus 

haut, associees aux modeles de localisation qui viennent 
d'etre' decrits, constituent, le modele polycristallin . Ce 
modele permet de calculer l'effet de sollicitations 
mecaniques sur un objet, mais pour cela il doit etre 
10 integre a un modele de calcul par elements finis. 

On va maintenant deer ire 1 1 integration du modele 
polycristallin de comporternent du materiau dans un modele 
de calcul par elements finis. 

On va tout d'abord decrire le cas d ! un calcul sans 
15 endommagement , e'est-a-dire dans lequel le modele- de 
comporternent ne comport e que des briques "grains" . Dans 
un tel modele, le nombre de briques "grains" etant 
faible, on affecte une seule brique "grain" a chaque 
point de discretisation du calcul de structure, par 
20 elements finis. Les briques "grains" sont affectees aux 
points de Gauss, qui sont les points d 1 integration dans 
le calcul par elements finis. En chacun -de ces points, le 
nombre ,de briques "grains" etant egal a un, seules les 
relations scalaires du comporternent des briques "grains" 
25 sont utilisees. De plus, les vitesses de deformation 
plastique macroscopiques sont confondues avec celles des 
briques "grain" . II en est de meme de la contrainte 
macroscopique qui est egale a la contrainte dans la 
brique "grain" et de la deformation macroscopique qui est 
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egale a la deformation dans la brique "grain" . II en 
resulte que les relations de localisation et 
d'homogeneisation du modele de changement d'echelle 
deviennent inutiles . 
5 Afin de bien representer la structure du materiau et 

en particulier son homogeneite dans l'espace, les briques 
"grain" sont affectees aux points d 1 integration du 
maillage de ' la structure au moyen d'un algorithme 
garantissant une permutation au sein d l un m§me element 
10 fini et entre les elements finis voisins. Ainsi, les 
briques "grain" sont reparties de fagon analogue aux 
grains du materiau reel . Un exemple de repartition est 
donne a la figure 6 correspondant au cas d'un materiau a 
structure cubique tel que decrit precedemment . Cette 
15 figure represente un Element fini rep£re g^neralement par 
10 comportant six faces en forme de quadrilatere . Les 
sommets et les milieux des aretes correspondent a des 
nceuds du modele reperes 11, ce qui fait vingt nceuds . Les 
points d' integration ou points de Gauss, reper6s 12, sont 
20 situes a 1 1 interieur de l 1 element fini. Les briques 
"grain" correspondant aux sept orientations sont 
reparties comme 1 1 indiquent les chiffres 1 a 7 associes 
aux points de Gauss. Dans le cas ou I'objet etudie est 
une plaque ou une coque on represente cet objet par 
25 plusieurs couches d'epaisseur identique avec trois points 
d 1 integration en milieu et en limite de chaque couche, ce 
qui conduit a un nombre de points d 1 integration egal a 2 
p + 1. Comme on utilise de trois a cinq couches le nombre 
de points d 1 integration est de sept, neuf ou onze . Ainsi, 
30 la plaque ou la coque comporte toutes les briques "grain" 
du modele dans son epaisseur, comme un materiau reel. Une 
permutation garantit une bonne representation du materiau 
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reel dans les directions tangentes a la plaque ou a la 
coque . Une permutation garantit egalement une bonne 
representation moyenne du mat^riau reel dans l'<§paisseur 
de la plaque ou de la coque. La repartition des briques 
"grain" est representee §. la figure 7. L 1 element fini 20, 
comporte trois couches 21, 22 et 23. Les points 
d' integration sont regroup£s en quatre ensembles 24, 25, 
26 et 2 7 comportant chacun sept briques "grain" dont les 
orientations sont reperees de 1 a 7 . Dans certains cas il 
peut etre interessant d'affecter plusieurs briques 
"grain" a un des points d 1 integration voire toutes les 
briques "grain" du modele a chaque point d 1 integration. 
Dans ce cas, .les calculs sont un peu plus longs. 

En particulier, dans le calcul par elements finis 
avec Evolution de texture utilisant les dix briques 
"grain" precedemment definies pour un matejriau 
initialement isotrope, on affecte a certains points 
d' integration deux briques "grains" correspondant ..• aux 
grains dissocies 1 et 8, 2, et 9, 3 et 10. Si I 1 on utilise 
les dix briques "grains" supplement aires obtenues en 
substituant les normales aux plans de glissement n s et 
les directions de glissement \. , on affecte ces briques 
aux mimes points d 1 integration que les briques la 10. Le 
nombre de comportements distincts des points 
d 1 integration reste ainsi limite £ sept, comme dans le 
cas d'un calcul sans evolution de texture pour un 
materiau isotrope . 

Lorsqu'on considere le cas dans lequel les joints 
peuvent §tre endommages, c'est-a-dire lorsqu'on utilise 
un modele comportant des briques "joint de grain", les 
fractions volumiques des briques "joint de grain" etant 
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tres faibles, pour obtenir une bonne representation du 
matgriau on associe une brique "grain" et une brique 
"joint de grain" a chaque point d 1 integration du modele 
de calcul par elements finis, et on utilise un algorithme 
5 de permutation de telle sorte que les couples brique 
"grain" /brique "joint de grain" ne soient pas toujours 
les memes . Dans ce cas, les relations de changement 
d'echelle restent necessaires mais sont tres simplif iees . 
Ces relations sont : 

10 Z « f q O* + f-jCF 

= y =g J =j 

f g + f j - 1 

15 a = S + 2 p. (l-p)[Tj(l ® 1 ) + cc l]fi(s P -5 p ) 

a = S + 2n (l-p) [ti (1 ® 1) + a l]f g U P -* p ) 
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Dans le cas oil l ! on veut modeliser 1 ' endommagement 
ductile, c'est-a-dire lorsqu*on utilise un modele 
comportant des briques "grain" et une brique "trou", on 
associe une brique "grain" et une brique "trou" a chacun 
des points d 1 integration. Dans ce cas #galement, les 
relations de localisation sont tres simplifiees et ces 
relations sont : 



• p 



(l-f) e 

.P 

f *n 



1-/ = 



S + Dipai exp 



per, 
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A-vec ce modele, il est alors possible de calculer 
1' evolution entre un instant 0 et un instant t, des 
proprietes du materiau dont est constitue l f objet pour 
lequel on fait 1 1 etude des contraintes et des 
def ormations §. 1 1 etat macroscopique . Cependant, pour 
faire un tel calcul, il est necessaire de pouvoir 
determiner les differents parametres qui interviennent 
dans le modele. 

Pour determiner ces parametres, une methode 
d 1 identification connue en elle-meme, selon laquelle on 
commence d'abord par .rassembler dans une base de donnee 
relative au materiau des resultats d 1 essais de traction, 
d 1 essais traction/ tors ion ou encore d ' autres essais 
mecaniques tels que des essais de fluage. Puis on simule 
par le calcul les differents essais de la base de donnee 
en utilisant des valeurs a priori des parametres. Puis on 
compare les - resultats de ces simulations aux essais 
proprement dits et, en fonction de 1 1 ecart entre les 
resultats experimentaux et les resultats calcules, on 
deduit de nouvelles valeurs pour les parametres. Enfin, 
on reproduit ce processus jusqu'a ce que les ecarts entre 
les resultats experimentaux et les resultats simules 
soient inf6rieurs a une valeur definie a 1 1 avance . 

A titre d^xemple, on a applique - le modele selon 
1 invention au cas d T un tube en Zircaloy recristallise 
soumis a un effort de traction biaxee. Dans le repere du 
Laboratoire (L X/ L 2 , L3) , on a utilise des coordonn^es 
. cylindriques . 

Lie modele a ete identifie sur la base de deux essais 
de traction dans les orientations axiales zz et 
tangentielle 66 du tube et d'un essai de torsion z6 du 
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tube. Le resultat de 1 ' identification est represents a la 
figure 9. Dans cette figure, les triangles reperes 31 
correspondent a des points experimentaux, la courbe 32 
represente les resultats tels qu'ils sont obtenus en 
appliquant les criteres de von Mises avec un modele 
macroscopique classique de comportement du materiau, la 
courbe 33 represente les resultats obtenus par un modele 
polycristallin selon 1 ' art anterieur, et la courbe 34 
represente les resultats obtenus avec le modele selon 
1" invention. Ces courbes montrent que le procede selon 
1* invention conduit a une tres bonne concordance des 
resultats experimentaux et des resultats calcules pour 
des tractions dans toutes les directions du plan zz-66, 
ce qui n'est pas le cas de 1 ' application simple des 
criteres de von Mises, ni le cas de 1 ■ utilisation d'un 
modele polycristallin selon 1 ' art anterieur. 
L' identification a demande 2 heures de calcul sur une 
station de calcul standard. L ' identification du modele 
polycristallin selon 1 1 art anterieur est decrite dans 
"Philippe Geyer, Comportement elasto-viscoplastique de 
tubes en Zircaloy-4 : approche experimental et 
modelisation micromecanique, these de l'Ecole des Mines 
de Paris, 9 decembre 1999". Cette identification a 
necessity plusieurs jours de calcul sur un ensemble de 5 
a 10 stations de calcul travaillant simultanement . Le 
resultat est moins bon qu'avec le modele selon 
1 1 invention non pas parce que le module polycristallin 
selon l'art anterieur est moins bon que le modele selon 
1' invention, mais parce qu'il conduit a des calculs tres 
longs qui n'ont pas permis d'optimiser 1 1 identification 
des parametres. Cet exemple montre 1 ' inter§t du procede 
selon 1 1 invention qui reduit de f agon considerable les 
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temps de calcul et ainsi permet de determiner I'effet de 
soil ici tat ions mecaniques sur des objets dans des 
conditions dans lesquelles les methodes selon 1 1 art 
anterieur ne sont pas applicables. 
5 Le procede qui vient d'etre decrit permet de calculer 

Involution des proprietes du materiau ainsi que ses 
endommagements . II fait appel i. des modeles visco- 
plastiques, mais I'homme du metier pourra utiliser tout 
modele de comportement anelastique compatible avec un 
10 modile polycristallin conforme a 1* invention et 
notamment, les comportements en fluage et en fatigue 
fluage. 

Les endommagements correspondent soit a 1 1 ouverture 
des briques 11 joint de grain 11 soit au grossissement, des 

15 briques "trou" . selon qu'on s » interesse a un endomraagement 
intergranulaire ou a un endommagement ductile. Dans;, les 
deux cas, l ! homme de metier peut definir un critere de 
rupture defini, par exemple, par un seuil d 1 endommagement 
tel, qu ! au dela, il y a rupture. On peut ainsi prevoir 

20 les conditions de rupture en service. Dans le cas de 
1 1 endommagement intergranulaire , dans le cadre d 1 une 
methode par elements finis, on peut definir un critere 
d f ouverture d'une brique "joint de grain" tel que, 
lorsque. ce critere est satisfait, il y a fissure et, au 

25 point correspondant, la resistance mecanique est nulle. 
Dans le cas de 1 1 endommagement ductile, dans le cadre 
d'une methode par elements finis, on peut definir une 
fraction volumique critique d'une brique "trou" telle 
que, lorsque ce critere est satisfait, il y a fissure et, 

30 au point correspondant, la resistance mecanique est 
nulle. Dans les calculs, on considere le point auquel ce 
critere est atteint en premier. Puis en prenant en compte 
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les points voisins sur lesquels les contraintes vont se 
reporter, on peut suivre la croissance de la fissure et 
ainsi determiner les conditions d' apparition d'une 
rupture . 

Le procede selon 1' invention peut §tre utilise pour 
determiner 1 1 incidence du procede de fabrication d'un 
objet comportant une operation de mise en forme par 
deformation plastique a froid, tels que 1 1 emboutissage, 
et en particulier pour determiner les proprietes 
mecaniques de 1' objet en chacun de ses points. En effet, 
1' incidence d'une operation de mise en forme par 
deformation plastique sur les proprietes mecaniques 
depend de la deformation ef f ectivement realisee et cette 
deformation depend en general du point considered Les 
proprietes mecaniques apres mise en forme ne sont done 
pas, en general, identiques en tous points de 1' objet. On 
peut alors utiliser devaluation de la repartition des 
proprietes mecaniques de 1" objet obtenue par le procede 
selon 1 • invention, pour determiner le comportement en 
20 service de cet objet. L'avantage du procede selon 
1' invention est que cette evaluation de la repartition 
des proprietes mecaniques d'un objet peut etre faite dans 
des conditions tres di verses a partir d'un nombre 
restreint d'essais mecaniques. En effet, pour evaluer les 
proprietes mecaniques d ' objet s constitues du meme 
mat£riau, mais realises par mise en forme selon des 
modalites diverses, il suffit de realiser les essais 
mecaniques necessaires a la construction du modele 
polycristallin selon 1' invention, representatif du 
materiau, puis appliquer le procede selon 1 ' invention a 
chacun des cas particuliers envisages. 



25 



30 



41 



A titre d'exemple, on peut utiliser le procede selon 
1 1 invention ■ pour evaluer le comportement au choc d'une 
piece de vehicules automobiles obtenue par emboutissage 
ou par pliage ou par tout autre procede de mise en forme 
induisant des deformations plastiques. Pour cela, on 
const ruit le moddle polycristallin represent at if du 
materiau dont est constitute la piece. On cale le module 
polycristallin a 1 1 aide des essais mecaniques de traction, 
ou tract ion/torsioii ou cisaillement necessaires . Puis, a 
l'aide d ! une methode de calcul telle qu'une methode de 
c.alcul par elements finis, incorporant le model e 
polycristallin * representatif du materiau, on simule 
1 1 operation de mise. en forme ' et on determine les 
proprietes mecaniques' de la piece en chacun de ses 
points. A partir de ce resultat, ' et en utilisant une. 
methode de calcul adaptee connue de'l'homme du metier, on 
determine, par exemple, les effets sur les pieces , ,d * un 
choc resultant d ! une collision. 

Liorsqu'on souhaite determiner, le comportement en 
service de differentes pieces realisees dans le m§rae 
materiau, on utilise le meme module polycristallin 
representatif du materiau pour effectuer les. calculs de 
simulation de comportement. Ainsi, il n'est pas utile de 
faire des essais mecaniques pour chaque pidce. 

Le procede selon l 1 invention peut egalement etre 
utilise pour determiner le comportement d'une piece 
constitute d'un materiau pouvant presenter une 
transformation de phrase, c r est-&-dire un changement de 
structure cristallographique sous I'effet des 
sollicitations thermomtcaniques qu'il subit dans les 
conditions de service ou de mise en forme considerees. 
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C'est le cas, par exemple, d'un alliage a memoire de 
forme isotrope qui possede deux phases, l'une stable a 
haute temperature (1 ' austenite) , 1' autre metas table a 
basse temperature (la martensite) . La diminution de la 
temperature et/ou 1 ' application d'une sollicitation 
mecanique induit un passage progressif de 1 1 austenite 
vers la martensite. La transformation inverse est 
possible de facon reversible. Cette transformation de 
phases s'accompagne d'une deformation anelastique qui 
s'ajoute a la deformation propre de chaque phase prise 
independamment . 

La simulation de cette transformation de phases est 
realisee en utilisant des briques "grains" (7 en 
1' occurrence pour un materiau isotrope), chacune etant 
orientee dans le repere local macroscopique (XI, X2, X3) 
par ses angles d'Euler W, cf> et <p . Le comportement de 
chaque brique est simule par les six systemes de 
glissement precedemment decrits et leurs lois 
d'ecrouissage connues par ailleurs. Les parametres de ces 
lois de comportement correspondent dans un premier temps 
a ceux de 1' austenite. 

En cours de chargement thermomecanique, en fonction 
d'un critere bien dSfini et connu, par exemple par 
"Leclercq S., Lexcellent C, A general macroscopic 
description of thermomechanical behaviour of shape memory 
alloys, J. Mech. Phys . Solids, Volume 44, pp 953-980, 
1996", l'austenite est susceptible de se transformer en 
martensite, a une vitesse donnee . A chaque etape du 
calcul, on verifie si ce critere est satisfait, et on 
attribue aux briques "grain" pour lesquelles le critere 
est satisfait, les caracteristiques mecaniques de la 
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martensite (evolution des param^tres des lois de 
comportement- des systemes de glissement de la martensite, 
connue par ailleurs) . La deformation de transformation de 
phase qui resulte de cette evolution, apparalt alors 
5 lorsqu'on applique les regies de transition d'echelle du 
moddle. 

Le modele qui vient d'etre decrit est essentiellement 
un modele numerique qui peut etre mis en oeuvre par tout 
moyen approprie, en particulier, par ordinateur en 
10 1 1 integrant dans un code de calcul par 616ments finis 
connu en lui -meme • 

Ce modele peut etre utilise de la meme fagon dans 
toute methode de calcul de deformation de structure et, 
notamment, dans les methodes de calculs aux differences 
15 finis ou dans des methodes qui ne font intervenir q'un 
point de l'objet etudie comme c'est le cas dans la 
• modelisation dii laminage de produits plats. L.'homme du 
metier sait faire les adaptations necessaires qui ne 
portent que sur 1 1 integration du modele dans la mdthode 
20 de calcul - 

On a decrit les exemples d 1 application ' & des 
materiaux ayant des structures cristallographiques 
cubiques ou hexagonales. Mais 1 ' invention ne se limite 
pas a ces exemples et I'homme du metier saura I'appliquer 
25 a d'autres structures cristallines . et 3. des textures 
di verses . Elle s 1 applique egalement a des structures 
mixtes constitutes de cristaux de natures differentes. 
Pour cela, il suffit de prendre en compte des briques 
"grain" dont une partie correspond S. une structure 
30 cristalline et une autre partie correspond a une autre 
structure cristalline . 
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En resume, le procede selon 1 ■ invention, applique a 
un objet const itue d'un materiau determine et soumis 
pendant une periode de temps determinee a des 
solicitations mecaniques et eventuellement thermiques 
5 determinees, consiste h : 

choisir un modele de comport ement du type 
polycristallin conforme a 1" invention, adapte au 
materiau, 

- ajuster les parametres du modele de comport ement a 
10 partir d'un nombre restreint d'essais mecaniques 

convenabl ement choisis , 

- identifier l'effet des sollicitations mecaniques et 
eventuellement thermiques auxquelles 1* objet est soumis 
pendant la periode de temps consideree, 

- et calculer I'effet des sollicitations en chaque 
point de 1' objet en utilisant une methode de calcul 
incorporant le modele de comportement du materiau. 

Le procede peut etre applique a tout type de materiau 
dont le comportement peut Stre represents par un modele 
du type polycristallin, et en particulier a tout materiau 
mStallique . 

Les sollicitations auxquelles sont soumises l'objet, 
qui peuvent St re des efforts ou des deplac ement s repartis 
dans l'espace et dans le temps et qui peuvent etre 
completees par des sollicitations thermiques, peuvent 
§tre determinees par tout moyen et en particulier par des 
enregistrements de mesures. 

Les methodes de calcul utilisees par le procede sont des 
methodes de calcul numerique mises en ceuvre par 
30 ordinateur a 1 1 aide de logiciels adaptes . 
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REVEND I CAT I ONS 

1. ProcSde pour determiner les effets d'une 
sollicitation mecanique sur un objet constitue d'un 
5 materiau solide selon lequel on calcule les deformations 
et les contraintes engendrees par la sollicitation 
mecanique en une plurality de points de 1' objet au moyen 
d'une methode de calcul numerique en mode anelastique 
dans laquelle le comportement du materiau solide est 
10 represents par un modele de comportement microscopique 
polycristallin utilisant une pluralite de briques 
« grain » dont les deformations sont determinees a partir 
d'une pluralite de syst£mes de glissement propres au 
materiau solide, la trace du tenseur des deformations 
15 microscopiques etant nulle, chaque brique "grain" ayant 
une orientation et une fraction volumique, caracterise : en 
ce que : 

- le nombre de briques « grains » est inferieur *6u 
egal £ dix, 

20 - le nombre de syst ernes de glissement est inferieur 

ou Sgal a six, 

en outre, le modele microscopique utilise eventuellement : 
jusqu'a sept briques « joint- de grain » * ayant 

chacune une orientation et une fraction volumique, et 
25 dont les deformations sont determinees a partir d'au 

moins un syst erne d'ouverture , la trace du tenseur des 

deformations plastiques microscopiques de chaque brique 

« joint de grain » etant positive, 

au plus une brique « trou » ayant une fraction 
30 volumique variable, dont les deformations plastiques sont 

purement volumiques et dont la trace du tenseur des 

deformations plastiques microscopiques est positive. 
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2. Procede selon la revendication 1, caracterise en 
ce que les six systemes de glissement associes a une 
brique « grain », d6finis dans un repere orthonorme 
associe a la brique « grain » par les normales aux plans 
de glissement et par les directions de glissement, sont 
tels que le double de la matrice d' orientation associee & 
chaque systeme de glissement ne soit constitue que de 0, 
1 ou - 1 . 

3 . Procede selon la revendication 1 ou la 
revendication 2, caracterise en ce que pour un materiau 
isotrope a structure cubique, les deformations sont 
determinees a partir de deux families de trois systemes 
de glissement, une premiere famille correspondant aux 
faces du cube, et la deuxieme famille correspondant aux 
plans a 45° de ces faces, et en ce que le nombre de 
briques « grains » est egal a sept, les briques 
« grains » etant obtenues par symetrisation de deux 
briques « grains » primitives situees sur les grands 
cercles a 45° de la figure de p61es. 

4. Procede selon l'une quelconque des revendications 
1 S. 3 caracterise en ce que on determine en outre 
Involution de la texture du materiau en calculant la 
rotation Q du reseau cristallin de chaque brique "grain" 

par rapport & un repere corotationnel . 

5 . Procede selon la revendication 4 caracterise en ce 
que, pour un materiau isotrope & structure cubique, le 
nombre de briques "grain" est <§gal a dix. 
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6. Proc<5de selon 1 1 une quelconque des revendication 1 
a 3, caracterise en ce que on determine en outre 
1 1 avancement d 1 une transformation de phase . 

7. Procede selon la revendication 6, caracterise en 
ce que on attribue a chaque brique "grain" , des 
caracteristiques fonction de 1 ' avancement de la 
transformation de phase. 



8. Procede selon I'une quelconque des revendi cat ions 
1 & 7, caracterise en ce que les orientations et les 
fractions volumiques des briques « grains » sont 
determinees par ajustement a partir d'essais de tractioft 
15 biaxiale et de traction torsion. 



9. Procede selon l'une quelconque des revendications 
1 & 8, caracterise en ce que, en 1 ' absence de brique. 
« trou », le modele de passage du niveau macroscopique au 
niveau microscopique pour une brique « grain », 
permettant de calculer le tenseur des contraiiites 
microscopiques a, dans la brique « grain » en fonction du 



tenseur des contraintes macroscopiques £ , du tenseur s P 

=g 

des deformations plastiques microscopiques dans la brique 
25 « grain » et du tenseur E P des deformations plastiques 
macroscopiques > s ' ecrit : 

£ s - £ + 2m(1 - /?) Ei? (1 ® 1) + a 1] (E P - s P ) 

avec, l/a=l + DJB~ , et 1/17 = 2 + 5*—, B~ est la valeur 

maximale atteinte par le deuxieme invariant de von Mises 
30 lors de la deformation. * 
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10. Procede selon l'une quelconque des revendications 
1 a 9, caracterise en ce que la methode de calcul 
numerique est une methode de calcul par elements finis et 
en ce que a chaque point d' integration de la methode de 
calcul par elements finis, on associe une seule brique 
« grain » et eventuellement une brique « joint de grain » 
ou une brique « trou », de telle sorte que la repartition 
des briques soit homogene. 

11. Proc€de selon l'une quelconque des revendications 
1 a 10, caracterise en ce que la methode de calcul 
numerique est une methode de calcul mecanique ou 
thermomecanique, faisant intervenir le temps. 

12. Procede selon la revendication 11, caracterise en 
ce que on utilise une brique « trou » et qu'a 1' instant 
initial du calcul numerique, la fraction volumique de la 
brique « trou » n'est pas nulle. 

13. Procede selon l'une quelconque des revendications 
1 a 11, caracterise en ce que le notnbre de briques "joint 
de grain" est superieur ou <§gal a un, et en ce que on 
definit un critere tel ' que, lorsqu'en un point de 
I'objet, le critere est satisfait, on attribue au point 
considere une resistance mecanique du materiau nulle, de 
fa?on a simuler la presence d l une fissure. 

14. Procede selon l'une quelconque des revendications 
1 a 8 et 10 a 12, caracterise en ce que on utilise une 
brique "trou", et en ce que on definit un critere tel 
que, lorsqu'en un point de I'objet, le critere est 
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satisfait, on attribue au point considere une resistance 
mecanique du mat<§riau nulle, de facon a simuler la 
presence d'une fissure. 

5 15. Procede selon l'une guelconque des revendications 

1 a 14, caracterisg en ce que l'objet en materiau solide 
est une piece metallique et en ce que la sollicitation 
mecanique est la sollicitation mecanique associee a une 
operation de mise en forme par deformation plastique 
10 telle ' que notamment 1' emboutissage, l'etirage, le 
laminage, le cintrage ou le forgeage. 

16. Procede selon l'une quelconque des revendications 
1 a 15, caracterise en ce que l'objet en materiau solide 

15 est. une piece metallique et en ce que la sollicitation 
mecanique est une sollicitation mecanique correspondant <a 
1' utilisation de la piece metallique. v 

17. Precede selon la revendication 16, caracterise en 
20 ce que 1 ' utilisation en service de la piece engendre au 

moins un phenomene, pris parmi les phenom^nes de fluage, 
de fatigue et de fatigue -fluage . 
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